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Abstract. We consider the problem of writing real polynomials as determinants of sym- 
metric linear matrix polynomials. This problem of algebraic geometry, whose roots go back 
to the nineteenth century, has recently received new attention from the viewpoint of convex 
optimization. We relate the question to sums of squares decompositions of a certain Her- 
mite matrix. If some power of a polynomial admits a definite determinantal representation, 
then its Hermite matrix is a sum of squares. Conversely, we show how a determinantal 
representation can sometimes be constructed from a sums-of-squares decomposition of the 
Hermite matrix. We finally show that definite determinantal representations always exist, 
if one allows for denominators. 



Introduction 

A polynomial p e M[x] in n variables x — (x±, . . . , x n ) with p(0) = 1 is called a real-zero 
polynomial if p has only real zeros along every line through the origin. The terms hyperbolic 
or real stable polynomial are also common and mean essentially the same, but usually for 
homogeneous polynomials. The typical example is a polynomial given by a definite (linear 
symmetric) determinantal representation 

p = det(I + A 1 x 1 H h A n x n ), 

where A\, . . . , A n are real symmetric matrices and I is the identity. A representation of this 
form is a certificate for being a real- zero polynomial. In other words, the fact that p is a real- 
zero polynomial is apparent from the representation. A definite determinantal representation 
also provides a description of the rigidly convex region of p. This is the closed connected 
component of the origin in the complement of the zero-set of p. It is always convex, and 
given a definite determinantal representation of p, it coincides with the set of points where 
the matrix polynomial / + A\Xi + • • • + A n x n is positive semidefinite. 
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In recent years, real-zero polynomials and their determinantal representations have been 
studied mostly with a view towards convex optimization, specifically semidefinite and hy- 
perbolic programming. In general, one would like to answer the following questions: 

(1) Under what conditions does a real- zero polynomial have a definite determinantal 
representation? 

(2) If such a representation exists, what is the minimal matrix dimension and how can 
the representation be computed effectively? 

(3) If no such representation exists, what other certificates for being a real-zero polyno- 
mial are available? 

Question (1) is the most immediate and has consequently received the most attention. It 
ties in with the theory of determinantal hypersurfaces in complex algebraic geometry, whose 
roots go back to the nineteenth century. Arguably the most important modern results are 
the Helton- Vinnikov theorem in [7], which gives a positive answer for n — 2, and Branden's 
negative results in higher dimensions in [4]. Since there are various subtle variations of the 
question, it is not always easy to figure out what is known and what is not; we give a very 
brief overview after the introduction below. 

Question (2), which should be of interest for practical purposes, has not been studied 
very systematically so far. Even in the case n = 2, the classical approach of Dixon for 
constructing determinantal representations is quite algorithmic in nature but hard to carry 
out in practice (see [5], or [14] for a more modern presentation.) 

One approach to Question (3) is to study the determinantal representability of a suitable 
power or multiple of p if no representation for p exists. This is motivated by the Generalized 
Lax Conjecture, as described below. On the other hand, the real-zero property does not 
have to be expressed by a determinantal representation. That a polynomial p in one variable 
has only real roots is equivalent to its Hermite matrix being positive semidefinite. This is 
a symmetric real matrix associated with p, which provides one of the classical methods for 
root counting. To treat the multivariate case, we use a parametrized version of the Hermite 
matrix with polynomial entries. In a typical sums-of-squares-relaxation approach common 
in polynomial optimization, we then ask for the parametrized Hermite matrix "H(p) to be a 
sum of squares, which means that there exists a matrix Q such that T-L(p) = Q T Q. (This is 
called a sum of squares rather than a square, because Q is allowed to be rectangular of any 
size). This approach has been used before by Henrion in [9] and by Parrilo (unpublished) as 
a relaxation for the real-zero property, which is exact in the two-dimensional case. 

While the Hermite matrix provides a practical way of certifying the real-zero property, 
having a definite determinantal representation of p is clearly much more desirable, since it 
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also yields a description of the rigidly convex region by a linear matrix inequality And even 
if one is only interested in the real-zero property, the multivariate Hermite matrix is a fairly 
unwieldy object compared to the original polynomial, and a sum-of-squares decomposition 
even more so. 

Our main goal is therefore to use a sum-of-squares decomposition of the parametrized 
Hermite matrix of a polynomial p to construct, as explicitly as possible, a definite determi- 
nantal representation of p, or at least of some multiple of p. We first show in Section 1 that 
a definite determinantal representation of some power of p of the correct size always yields 
a sum-of-squares decomposition of H(p) (Thm. 1.6). In Section 2, we make an attempt at 
the converse. This is partly motivated by our experimental finding that the Hermite matrix 
of the Vamos polynomial, which is the counterexample of Branden, is not a sum of squares 
(Example 1.9). Note also that in the case n = 2, where every real- zero polynomial possesses 
a definite determinantal representation by the Helton- Vinnikov theorem, the parametrized 
Hermite matrix can be reduced to the univariate case. It is therefore a sum of squares if 
and only if it is positive semidefinite, by a result of Jakubovic [10]. Given a decomposition 
'Hip) = Q T Q, we show that a definite determinantal representation of a multiple of p can 
be found if a certain extension problem for linear maps on free graded modules derived 
from Q has a solution (Thm. 2.5). Given Q, the search for such a solution amounts only 
to solving a system of linear equations. This method can in principle also be applied if 
the sums of squares decomposition uses denominators. Finally, we show that by allowing 
a sum-of-squares decomposition with denominators, which exists whenever H{p) is positive 
semidefinite, one can always obtain a determinantal representation with denominators: 

Theorem. Let p be a square-free real-zero polynomial with p(0) = 1. There exists a sym- 
metric matrix M. whose entries are real homogeneous rational functions of degree 1 such that 
p = det(I + M). 

The precise statement is given in Thm. 3.1. 

Acknowledgements. We would like to thank Didier Henrion, Pablo Parrilo, Rainer Sinn, 
and Cynthia Vinzant for helpful comments and discussions. 

Known results 

o For n — 2, every real-zero polynomial of degree d has a real definite determinantal 
representation of matrix size d by the Helton- Vinnikov theorem [7]. 

o For n > 3 and d sufficiently large, a simple count of parameters shows that only 
an exceptional set of polynomials can have a real determinantal representation of 
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size d. The question whether every real-zero polynomial has a definite determinantal 
representation of any size became known as the generalized Lax conjecture. 

o The generalized Lax conjecture was disproven by Branden who even showed the exis- 
tence of real-zero polynomials p such that no power p r has a determinantal represen- 
tation of any size [4]. His smallest counterexample, the so-called Vdmos polynomial, 
is of degree 4 in 8 variables (see 1.9 below). 

o Netzer and Thorn [11] have proved that only an exceptional set of polynomials can 
have a determinantal representation, even if one allows for matrices of arbitrary size. 
This is true for n > 3 and d sufficiently large, or d > 4 and n sufficiently large. They 
also show that if p is a real-zero polynomial of degree 2, then there exists r > 1 such 
that p r has a determinantal representation. On the other hand, there exists such p 
where one cannot take r = 1. 

o Another result of Helton, McCullough and Vinnikov [8] (see also Quarez [13]) says 
that every real polynomial has a real symmetric determinantal representation, though 
not necessarily a definite one. This means that the constant term in the matrix 
polynomial cannot be chosen to be the identity matrix in their result. 

o The most general form of the Lax conjecture says that every rigidly convex set is 
a spectrahedron. In terms of determinantal representations, this amounts to the 
following: For every real- zero polynomial p there exists another real-zero polynomial 
q such pq has a real definite determinantal representation and such that q is non- 
negative on the rigidly convex set of p. This conjecture is still wide open, even without 
the additional positivity condition on q. Note that if pq has a definite determinantal 
representation, then q is automatically a real-zero polynomial. 



1. The Hermite Matrix 

In this section we introduce the parametrized Hermite matrix H(p) of a polynomial. It 
is positive semidefinite at each point if and only if p is a real zero polynomial. If some power 
of p admits a determinantal representation of the correct size, then H(p) even turns out to 
be a sum of squares of polynomial matrices. 

Let p = t d + p x t d ~ x + ■■■+ p d ^t + p d e R[t] be a monic univariate polynomial of degree 
d and let A i, ... A^ be the complex zeros of p. Then 

d 

N k (p) = J2^ 

i=l 
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is called the k-th Newton sum of p. The Newton sums are symmetric functions in the roots, 
and can thus be expressed as polynomials in the coefficients p^ of p. The Hermite matrix of 
p is the symmetric d x d matrix 

H(p) : (.V;., ,(/>))„, , r 

It is a Hankel matrix whose entries are polynomial expressions in the coefficients of p. Note 
that Hip) = V T V, where V is the Vandermonde matrix with coefficients Ai, . . . , \d- 

The following well-known fact goes back to Hermite. For a proof, see for example 
Theorem 4.59 in Basu, Pollack and Roy [1]. 

Theorem 1.1. Let p G R[i] be a monic polynomial. The rank of H(p) is equal to the number 
of distinct zeros of p in C. The signature of the Hermite matrix H(p) is equal to the number 
of distinct real zeros of p. 

In particular, H(p) is positive definite if and only if all zeros of p are real and distinct, 
and H(p) is positive semidefinite if an only of all zeros are real. □ 

Now let p G be a polynomial of degree d in n variables x = (xi,. . . ,x n ). The 
polynomial p is called a real-zero polynomial (with respect to the origin) if p(0) = 1 and 
for every a G M. n , the univariate polynomial p(ta) G has only real zeros. We want to 
express this condition in terms of a Hermite matrix. Write p = Ylt=o Pi with Pi homogeneous 
of degree i, and let P(x,t) = Yli=oPi^ d ~ l be the homogenization of p with respect to an 
additional variable t. We consider P as a monic univariate polynomial in t and call the 
Hermite matrix H(P) the parametrized Hermite matrix of p, denoted H(p). Its entries are 
polynomials in the homogeneous parts pi of p. The (i, j)-entry is a homogeneous polynomial 
in x of degree i + j — 2. 

Corollary 1.2. A polynomial p G M.[x] with p(0) = 1 is a real-zero polynomial if and only if 
the matrix H(p)(a) is positive semidefinite for all a6l™. 

Proof. By Theorem 1.1, H(p)(a) is positive semidefinite for a G M n if and only if the uni- 
variate polynomial t d p(ait _1 , . . . , a n t— 1) has only real zeros. Substituting t~ l for t, we see 
that this is equivalent to p(ta) having only real zeros. □ 

The following is Proposition 2.1 in Netzer and Thorn [11]. 

Proposition 1.3. Let M. = x\Mi + - ■ -+x n M n be a symmetric linear matrix polynomial, and 
let p = det(I — Ai). Then for each a G R n , the nonzero eigenvalues of M.{a) are in one to 
one correspondence with the zeros of the univariate polynomial p(ta) , counting multiplicities. 
The correspondence is given by the rule A t— > y. □ 
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Lemma 1.4. Let p £ M[x] be a real-zero polynomial of degree d, and assume that p r = 
det(I — M) is a symmetric determinantal representation of size k, for some r > 0. Then 

^W = ^(tr(^- 2 )), 
except possibly for = (1, 1), where H(p)i,i = d and tr(.M ) = k. 

Proof. For each a £ M. n , the trace of Ai(a) s is the s-power sum of the nonzero eigenvalues 
of Ai(a). These eigenvalues are the inverses of the zeros of p(ta), by Proposition 1.3, but 
each such zero gives rise to r many eigenvalues. Since the zeros of pita) correspond to the 
inverses of the zeros of t d p(t~ 1 a), the trace of A4(a) s equals the s-power sum of the zeros of 
t d p(t~ 1 a) multiplied with r. This proves the claim. □ 

Definition 1.5. Let H £ Sym d (R[x]) be a symmetric matrix with polynomial entries. H is 
a sum of squares, if there is a d' x d- matrix Q with polynomial entries, such that 7i = Q T Q. 
This is equivalent to the existence of d' many d- vectors Qi with polynomial entries, such 
thatft = £Li QiQj. 

Theorem 1.6. Let p £ M[x] be a real-zero polynomial of degree d. If a power p r admits 
a definite determinantal representation of size r ■ d, for some r > 0, then the parametrized 
Hermite matrix Hip) is a sum of squares. 

Proof. Let p r = det(I — M) with Ai of size k = rd, and denote by the (£, m)-entry of 
M s . Put Q im = (qj®, q^Y £ R[x] d . Then we find 



k 



E G*»QL = E = ("•(-^ ; M' ')), ; , = rH( P ), 

e ' m=1 W=i / ,. ; i (/ 

by Lemma 1.4. □ 

Remarks 1.7. (1) If the determinantal representation of p r is of size k > rd, then H(p) 
becomes a sum of squares after increasing the (1, 1) -entry from d to k/r. This is clear from 
the above proof. 

(2) It was shown in Netzer and Thorn [11] that if a polynomial p admits a definite 
determinantal representation, then it admits one of size dn, where d is the degree of p and 
n is the number of variables. So if any power p r admits a determinantal representation of 
any size, then ~H(p) is a sum of squares, after increasing the (1, l)-entry from d to dn. Note 
that this is independent of r. 

(3) The determinant of H(p) is the discriminant of t d p(t~ 1 x) in t. If H(p) = Q T Q, it 
follows from the Cauchy-Binet formula that the determinant of H(p) is a sum of squares in 
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R[x]. Thus, by the above theorem, the discriminant of det(tJ + M.) in t is a sum of squares, 
a fact that has long been known, at least since Borchardt's work from 1846 [3]. 

(4) The sums-of-squares decomposition of Hip) obtained by Thm. 1.6 from a deter- 
minantal representation p r = det(J — Ai) is extremely special. In principle, it is possible 
to characterize the decompositions of H{jp) coming from a determinantal representation by 
a recurrence relation that they must satisfy. But this does not appear to be a promising 
approach for finding determinantal representations. 

Example 1.8. It was shown in Netzer and Thorn [11] that if p is quadratic, a high enough 
power admits a definite determinantal representation of the correct size. Thus Hip) is a sum 
of squares in this case. This can also be shown directly. Write 



P 



x T Ax + b T x + 1 



with A e Sym n (IR) and b € M. n . Then p is a real- zero polynomial if and only if bb T — AA >z 0, 
as is easily checked. We find t 2 p(t~ 1 x) = x T Ax + b T x -t + t 2 , and so we compute 

2 -b T x 
-b T x x T (bb T -2A)x 

Write bb T - AA 



Hip) 



EI 



=1 ViV- 



as a sum of squares of column vectors Vi G 



Set 



/ 1 -hb T x \ 



Q 



1 T 
2^1 X 



\v T n x ) 



Then Uip) = 2 • Q T Q. 



□ 



Example 1.9. We consider Branden's example from [4]. It is constructed from the Vamos 
cube as shown in Figure 1. Its set of bases B consists of all four element subsets of {1, . . . , 8} 
that do not lie in one of the five affine hyperplanes. Define 

a degree four polynomial in IR[xi, . . . , £§]. It contains as its terms the product of any choice 
of four pairwisely different variables, except for the following five: 



X1X4X5XQ, X2X3X5XQ, X2X^XjXg, £1X4X7X8, X1X2X3X4. 

Now p = q(xi + 1, . . . , Xg + 1) turns out to be a real-zero polynomial, of which Branden has 
shown that no power has a determinantal representation. 
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Figure 1. The Vamos Cube 

6 



1 




3 



7 



We can apply the sums-of-squares-test to the Hermite matrix "H(p) here. Unfortunately, 
the matrix is too complicated to do the computations by hand. When using a numerical 
sums-of-squares-plugin for matlab, such as Yalmip, the result however indicates that T-L(p) is 
not a sum of squares. In view of Theorem 1.6 this shows again that no power of p admits a 
determinantal representation. Note that if some power p r has a determinantal representation, 
then it has one of size 4r. This was proven by Branden or follows more generally from Netzer 
and Thom, Theorem 2.7 [11]. 

Finally, we can apply the sums-of-squares-test also to small perturbations of Branden's 
polynomial. For example, p can be approximated as closely as desired by real-zero poly- 
nomials, which have only simple roots on each line through the origin (in other words, the 
Hermite matrix is positive definite at each point o ^ 0). Such a smoothening procedure is 
for example describe in Nuij [12]. Still, Yalmip reports that the Hermite matrix is not a sum 
of squares, if the approximation is close enough. This is exactly what one expects, since the 
cone of sums of squares of polynomial matrices is closed, and the Hermite matrix depends 
continuously on the polynomial. 



In this section we are interested in the converse of the above result. Namely, can a sums- 
of-squares decomposition of T-L(p) be used to produce a definite determinantal representation 
of p or some multiple? We describe a method to do this, which amounts to only solving a 
system of linear equations. 

Let p = 1 +pi + ■ ■ ■ +Pd £ $L[x] be a real-zero polynomial of degree d. Since the matrix 
%(jp) is everywhere positive semidefinite, it can be expressed as a sum of squares if one allows 
denominators in M.[x\. This generalization of Artin's solution to Hilbert's 17th problem was 
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first proved by Gondard and Ribenboim in [6]. We need to make a slight adjustment to our 
situation. 

Lemma 2.1. There exist a matrix polynomial Q G Matfc X d(lR[x]) , for some k > 0, and a 
homogeneous non-zero polynomial q G M[x] such that 

q 2 U(p) = Q T Q. 

Proof. By the original result of Gondard and Ribenboim [6] there is some non-zero polyno- 
mial q G such that q 2, H{p) = Q T Q for some Q G Matfc X( f(R[a;]). We want to make q 
homogeneous. 

Write q = q r + q r+ \ + ■ ■ ■ + qa, where each qi is homogeneous of degree i, and q r ^ 0, 
qR 7^ 0. Since the i-th diagonal entry in 1-L{p) is homogeneous of degree 2(i — 1), each entry 
in the i-th column of Q has homogeneous parts of degree between r + i — 1 and R + i — 1. Let 
Qmin be the matrix one obtains from Q by choosing only the homogeneous part of degree 
r + i — 1 of each entry in each i-th column. Put Q = Q — Q m i Q and note that all entries in the 
i-th column of Q have non-zero homogeneous parts only in degrees at least r + i. We now 
compute q 2 H(p) = Qj nin Q m m + Qmin Q + 2 T Qmin + Q T Q, compare degrees on both sides, 
and find qll-L(p) = <2 min Q m in, as desired. □ 

We will now describe the setup that we are going to use for the rest of this section. We 
fix a representation of q 2 l-i{p) = Q T Q as in Lemma 2.1. As before, let P = t d ■ p(t~ 1 x) = 
t d + pit 11 ' 1 + • • • + pa G R[x, t], and consider the free ]R[x]-module 

d-l 

A = R[x, t]/(P) = R[x] ■ f ^ R[x] d . 

i=0 

Since P is homogeneous, the standard grading induces a grading on A. We shift this grading 
by r, the degree of q, and obtain a grading with deg(t l ) = r + i for i = 0, . ..,d — 1. 
This turns A into a graded K [x] -module, where WL[x] is equipped with the standard grading. 
Furthermore, we equip A with a symmetric M.[x] -bilinear and IR[a;]-valued map (-, -) p defined 
by 

(f,g) p --=f T {q 2 n P ))g, 
for / = (f 1} . . . , f d ) T and g = {g u . . . , g d ) T in A. 

Next, consider the map Ct : A — > A given by multiplication with t. This is an -linear 
map which we can compute with respect to our chosen basis: 

A : (A, • • • , fd) T H> (-Pdfd, fi - Pd-ifd, fd-l - Pifd) T - 
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Note that C t is of degree 1 with respect to the grading, i.e. deg(£ t (/)) = deg(/) + I. We 
identify C t with the matrix that represents it, so that 

-Pd ^ 

-Pd-i 

1 

-Pi J 

which is exactly the companion matrix of P, viewed as a univariate polynomial in t. It is 
well known and easy to see that P is the characteristic polynomial of C t , so that 

det (/ -C t )=p 

Lemma 2.2. The linear map C t is self-adjoint with respect to (•, -) p , i.e. 

(£tf,g) P = (f,£ t g) P 

holds for all f,g e A. 

Proof. We may divide by q 2 on both sides and hence assume that q — 1. It is enough to 
show (C t ei,ej) p = (ei,£ t Zj)p f° r all i,j, where C{ is the i-th unit vector. For i,j < d, this 
follows from the fact that 7i(p) is a Hankel matrix. For % — j — d, it is clear from symmetry. 
So assume j < i = d. We find 

d d 
(£ t e d ,ej) p = -^p d ^ i+1 e i n(p)e j = -^2p d _ i+1 N i+j _ 2 , 
i=i i=i 

where Nk is the k-th. Newton sum of P. On the other hand, we compute (e d , C-t£j) P = 
(e d , e.j+i) P = N d+ j-i. In conclusion, we have to show that 

d 

^pd-iNi+j^ = 0, 

i=0 

where we have set p — 1. This statement is equivalent to "Ylt^Pi^k-i — 0, where k = 
d + j — 1 > d. This last equation, however, follows immediately from the Newton identity 
kpk + YliZo Pi^k-i — 0, where we let p k = for k > d. □ 

Let B = M\x\ k . The k x ci-matrix Q in the decomposition of H(p) describes an R[x]- 
linear map A = R[x] d — > B, / i->- Qf. From the degree structure of V.(p), we see that each 
entry in the i-th column of Q is homogeneous of degree r + i — 1. So Q is of degree with 
respect to the canonical grading on B. 



( 

1 

o '•• 

v° ••• 1 



Lemma 2.3. 
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(1) If p is square-free, then Q: A — >■ B is infective. 

(2) We have 

(f,g) P = (Qf, Qg) 

for all f,g G A. In other words, Q is an isometry, taking ( , ) p to the canonical 
bilinear form ( , } on B. 

Proof. (2) is immediate from the fact that q 2 T-i{p) = Q T Q- (1) If Qf = 0, then 

0=(Qf, Qf) = (fJ) p = q*.f T U(p)f. 

For each a G W l for which p(ta) has only distinct roots, the matrix H(p)(a) is positive 
definite. So f(a) = for generic a, and thus / = 0. □ 

Time for a brief summary of what we have done so far. 

Setup 2.4. 

o Let p G R[x] be a real-zero polynomial of degree d with p(0) = 1, and let H(p) be 
its parametrized Hermite matrix. Fix a decomposition q 2 I-i{p) = Q T Q, where q is 
homogeneous of degree r and Q is a matrix of size k x d with entries in M[x\. 

o We have equipped the free module A = M.[x] d with a particular grading and with a 
bilinear form ( , ) p : A x A — > A. 

o Let B = M.[x] k be equipped with the canonical bilinear form and the canonical grad- 
ing. 

o The map Q : A — > B is an isometry and of degree 0. 

o Let C t be the companion matrix of t d p(t~ 1 x) with respect to t, so that 

det(7 - Ct) = p. 

The map C t : A — > A is self-adjoint with respect to (•, -) p and of degree 1. 
The following is our main result. 

Theorem 2.5. Let p G M[x] be a square-free real-zero polynomial of degree d with p(0) = 1. 
Assume that there exists a homogeneous symmetric linear matrix polynomial M. of size kx k 
such that the following diagram commutes: 

R[x] d = A — ^ B = R[x] k 



M 



R[x] d = A > B = R[x] k 

Then p divides det(I — M). 
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Remark 2.6. Note that the above described setup exactly means that we can hope for such 
a linear symmetric Ai to exist. Indeed the "strange" symmetry of £ t is transformed into the 
standard symmetry by Q, and the "strange" grading is translated to the standard grading. 

Proof. For generic a G M. n , the map Q(a) is injective by Lemma 2.3. Therefore, all eigenvalues 
of C t (a) are also eigenvalues of Ai(a). The eigenvalues of C t (a) are precisely the zeros of 
P(t, a), i.e. the inverses of the zeros of p(ta). So q = det(I — Ai) vanishes on the zero set of 
p, by Proposition 1.3. Since p is a square-free real zero polynomial, the ideal (p) generated 
by p in M[x] is real-radical (see Bochnak, Coste and Roy [2], Theorem 4.5.1(v)). It follows 
that q is contained in (p), in other words p divides q. □ 

Remark 2.7. Whether there exists such Ai can be decided by solving a system of linear 
equations. Indeed, set Ai = X\M\ + • • • + x n M n , where the Mj are symmetric matrices with 
indeterminate entries. The equation AiQ = QC t of matrix polynomials can be considered 
entrywise, and comparison of the coefficients in x gives rise to a system of linear equations 
in the entries of the M~. 



Example 2.8. Let p G M.[x] be quadratic. Write p = x T Ax + b T x + 1 with A G Sym n (IR) and 
b G M n . We have seen in Example 1.8 that H(p) admits a sums of squares decomposition if 
p is a real-zero polynomial, given by the matrix 

/ 1 





Q = V2 



.IhT 

2 
1, 



b T X \ 



1 T 
if i X 



1 T 
2 V n X 



J 



if bb T -4A = j:i 



It is now easy to find a homogeneous linear matrix polynomial Ai 



that makes the diagram in Theorem 2.5 commute, namely we can take 

M 

The resulting determinantal representation is 

det (J - M) 





/ 


-b T x 


T 
v{ X 




T 


\ 


1 




T 

v{ X 


-b T x 










2 ' 


















V 


T 








-b T x 


J 



H b T x 



n-l 



To give an explicit example, consider p = (xi + y/2)' 



x 2 



— Xa — xf, which itself 



does not admit a determinantal representation (by Netzer and Thorn [11]). The procedure 
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just described now gives rise to the linear matrix polynomial 



M = 



/ 


-V2X! 


Xl 


X2 


X3 


X4 


x 5 






Xl 


-V2 Xl 


















x 2 





-y/2xi 















x 3 








-y/2x! 












X4 











-V2xt 







\ 


x 5 














-V2xi 


J 



and finally 

det(I-M) = (l + V^x^-p. 

Example 2.9. There are also examples where no suitable M exists. We are grateful to 
Rainer Sinn and Cynthia Vinzant for helping us find this example. Consider the plane cubic 
p = (xi — l) 2 (xi + 1) — x\. One computes 

(3 xi 3x1 + 2x\ \ 

xi 3xl + 2x 2 2 x\ + 3xix 2 2 = Q T Q, 

3x\ + 2x\ x\ + 3xix% 3x\ + %xlx 2 2 + 2x\ ) 

where 

/ x 2 axix 2 \ 

n _ — x 2 6x1X2 

V2 V2xi V2(xl + x 2 2 ) 
\ 1 -xi x\ J 

and a = \{\fi + 1), b = \{\/7 — 1)- The equation MQ = QC t has 12 entries, each of 
which gives rise to several linear equations by comparing coefficients in x. One can check 
that already the equations obtained from the first two rows of MQ = QC t are unsolvable. 

3. RATIONAL REPRESENTATIONS OF DEGREE ONE 

There is always a way to make the diagram from the last section commute, if one allows 
for rational linear matrix polynomials. This will lead to rational determinantal representa- 
tions, as described now. 

Let p be a square-free real-zero polynomial. Since the parametrized Hermite matrix 
T-L{p) evaluated at a point a G lR n is positive definite for generic a, the matrix polynomial 
V.(p) is invertible over the function field R(x). Recall that the degree of a rational function 
/ /g E R(x) is defined as deg(/) — deg(g). Furthermore, we say that f/g is homogeneous 
if both / and g are homogeneous, not necessarily of the same degree. Equivalently, f/g is 
homogeneous of degree d if and only if (f/g)(Xa) = \ d (f/g)(a) holds for all a G lR n with 
9(a) ± 0. 
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Theorem 3.1. Let p be a square-free real-zero polynomial. Write q 2 l-i{p) = Q T Q with q 
homogeneous as in Lemma 2.1 and let 

M :=q- 2 QC t n(p)- 1 Q T . 

The matrix M. is symmetric with entries in M(x) homogeneous of degree 1, and satisfies 

det(7 - M) = p. 

Proof. Abbreviate TL(p) by "H and £ t by C. By Sylvesters determinant theorem, we have 
det(ifc — AB) = det(Id — BA) for any matrix polynomials A of size k x d and B of size d x k. 
In our situation, this yields 

det (4 - M) = det(/ fc - q- 2 QCV.- 1 Q T ) = det(I d - q~ 2 CHT X Q T Q) = det(I d - C) = p. 

We find 

M T = q- 2 Q{H- 1 ) T C T Q T = q- 2 QC'H- 1 Q T = M, 

where we have used C T % = % J £, which is Lemma 2.2. Thus M. is symmetric. 
Let r be the degree of q. By examining the degree structure of q 2 H, we find 

Q(\a) = Q(a) ■ diag(A r , X r+1 , . . . , A r+d_1 ) 

n(Xa) = diag(A°, . . . , A^ 1 ) • H{a) ■ diag(A°, . . . , A* -1 ) 

£{\a) = diag(A d , . . . , A 1 ) ■ £{a) ■ diag(A- rf+1 , X~ d+2 . . . , A ) 

for all a G lR n and A / 0. Hence for all a G lR n for which %{a) is invertible and q(a) ^ 0, 
and all A 7^ 0, we have 

M{Xa) = X- 2r q(a)- 2 Q(a) ■ diag(A r , . . . , A^" 1 ) • diag(A d , . . . , X)£{a) ■ diag(A- d+1 , . . . , A ) 
• diag(A°, . . . , X~ d+1 ) ■ H-\a) ■ diag(A°, . . . , X- d+1 ) ■ diag(A r , . . . , A^ 1 ) • Q T (a) 
= X- 2r q(a)- 2 Q(a)X r+d C(a)X- d+1 n- 1 (a)X r Q T (a) 
= X-M{a). 

□ 

Remark 3.2. Note that a representation p = det (J — Ai) as in Theorem 3.1 gives an alge- 
braic certificate for p being a real-zero polynomial. Since pita) = det (J — tAi(a)) , using 
homogeneity, the zeros of p(ta) are just the inverses of the eigenvalues of A4(a). Since M. 
is symmetric, all of these zeros are real. Theorem 3.1 now states that such an algebraic 
certificate exists for each real-zero polynomial p. 
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Example 3.3. Consider the quadratic polynomial p 



(x 1 + l) 2 — x\ — x\ — x\. We have 



n 



-2x 1 2{x\ 



-2x x 



Q T Q with Q 1 



V2 
— \f2xx y/2x 2 y/2xz \f2~X4 



which results in 



M 



f-Xi 

■1*2 

\ X4 



X2 ^ 

.(, .-, ~\~ X . j ~\~ X ^ 
X\ X2X4 



X3 

X1X2X3 



.7^ 



3-1 ^4 



\ 
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